An effective gluon potential and hybrid approach to Yang-Mills thermodynamics 
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We derive the partition function for tire SU(3) Yang-Mills theory in the presence of a uniform 
gluon field within the background field method. We show, that the n-body gluon contributions in 
the partition function are characterized solely by the Polyakov loop. We express the effective action 
through characters of different representations of the color gauge group resulting in a form deduced 
in the strong-coupling expansion. A striking feature of this potential is that at low temperature 
gluons are physically disfavored and therefore they do not yield the correct thermodynamics. We 
suggest a hybrid approach to Yang-Mills thermodynamics, combining the effective gluon potential 
with glueballs implemented as dilaton fields. 

PACS numbers: 12.38. Aw, 25.75. Nq, ll.lO.Wx 



1. INTRODUCTION 

The SU{Nc) pure gauge theory has a global Z{Nc) 
symmetry which is dynamically broken in the high tem- 
perature phase. The Polyakov loop, defined from the 
temporal gauge field integrated over the Euclidean time, 
plays a role of an order parameter of the Z {Nc)sloha\ 
symmetry [l| . Effective Polyakov- loop models 0, U have 
been suggested as a macroscopic approach to the pure 
SU(3) gauge theory. The thermodynamics that comes 
out from such models is qualitatively in agreement with 
that obtained in lattice gauge theories Q. Alternative 
approaches are based on the quasi-particle picture of 
thermal gluons A natural extension was carried out 
where the gluons propagate in a constant gluon back- 
ground and the Polyakov loop appears in the partition 
function reflecting the group characters [6l-[l0j. 

In previous studies the effective potentials for gluons 
were approximated into the Boltzmann distribution and 
expanded in series of the Polyakov loop to finite order. 
Such expansion is however unclear to be justified in the 
SU(3) gauge theory around the first-order phase transi- 
tion. 

In this work, our main scope is to derive the SU(3) 
gluon thermodynamics utilizing the background field 
method for the one-loop quantization and to clarify 
the Polyakov-loop effective potentials used in the liter- 
ature anchored to the field theoretical basis. We show, 
that in this approach the calculated gluon potential, ex- 
hibits the correct asymptotic behavior at high temper- 
atures, whereas at low temperatures, it disfavors gluons 
as appropriate dynamical degrees of freedom. We derive 
its correspondence to the strong-coupling expansion, of 
which the relevant coefficients of the gluon energy dis- 
tribution are specified solely by characters of the SU(3) 
group. 

The paper is organized as follows: In Section [2] we 
derive the thermodynamic potential for the SU(3) pure 
Yang-Mills theory in terms of the Polyakov loop in the 
fundamental representation. A correspondence to other 
phenomenological potentials is deduced in Section |31 In 



Section H] we introduce a hybrid approach that matches 
gluons with glueballs at deconfinement transition and 
study its thermodynamics. Our concluding remarks are 
given in Section [5l 



2. MODELING GLUONS IN HOT MATTER 

To formulate thermodynamics of the SU(3) Yang-Mills 
theory, we start with the partition function for gluon A^^ 
and ghost C fields. 



Z = y VAf^VCVC exp 



(2.1) 



with the gauge fixing (GF) and the Faddeev-Popov ghost 
(FP) terms. The kinetic term is given by 

^fiu ^ dfj^Aiy — dj^Afj^ — i [Afj^, Ajj] , (2-2) 

where = A°T^ and tr [TT^] = 5"V2. 

Following [2l.1T]| we employ the background field gauge 
to evaluate the functional integrals. The gluon field is de- 
composed into background (classical) and the quan- 
tum Afj^ contribution 



Af_i — Af^ 



gA^. (2.3) 
We fix the background field gauge with the following form 
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where a is the gauge fixing parameter and D^A'^ = 
i [Af^,A'^'\. Expanding the Lagrangian and col- 



d^^A'^ 



lecting the terms quadratic in the quantum fields one 
gets 



£(2) 
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where we define 



transformation, one finds 



all 
b 

fi,ab 
V"/3 



(D )^^^^ 



-2iiv[A^ [T„,^fc]]5"^ 
-4ztr[A"'^[r,,Tb]] , 
5abd^, + 2itr [A^ [Ta,Tb]] 



(2.6) 



Here 't Hooft-Fcynman gauge (a = 1) was taken . 
Note that the ghost term does not contain A^^ and there- 
fore the Gaussian integral over the ghost fields can easily 
be carried out. In the following we keep only the terms 
quadratic in and rewrite Eq. p.Sp as 



faccfcbd9°'^A'i^A^''^ + 2fabcA" 



I3,c 



Ap . 



(2.7) 



In the above Lagrangian we consider a constant uniform 
background Aq: 



A^ = A'^^S. 



(2.8) 



It is convenient to take a diagonal and traceless genera- 
tors, i.e. a = 3,8. We first consider the simplest case 
where only contributes. Then, in the momentum 
space, the Lagrangian (|2.7p is of the following form 



£(2) 



-Al (D-^) , 

2 " \ ' ab 



{D 1)^^ = SabP^ + 2ifab3Alpo - facafcbS {M) . 

(2.9) 

where non-vanishing elements of the inverse propagator 
are given by 



D- 



'33 



= D- 



'88 



55 



77 
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-{D-^)^^^2zAIp,, 

iAlpo . 



67 



(2.10) 



Diagonalizing D ^ into D ^ — WD using unitary 



The partition function must be independent of gauge, i.e. 
din Z/ da = 0. Since the running coupling depends also on a, 
the condition reads, dlnZ/da = {d/da + {dg / da){d / dg)) \nZ = 
I12II . In this paper we work in Feynman gauge since the stan- 
dard partition function of a free boson gas is readily obtained in 
the high temperature limit. 
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bf 



(2.11) 

The two elements, {a,b) — (3,8), are never mixed with 
the background field Aq because the generators and 
are commuting, [T^,T^] = 0. 
In Euclidean space we replace po and Aq with ia;„ = 
i2m:T and — z^4, respectively. Then, with the prop- 
agator (|2.1ip . the summation over n which appears in 
the partition function is easily performed and one arrives 
at #2 



-\p\/T 



, (2.12) 



^Indet {D~^) = Indet (l - L^f 

n 

with the matrix, 

i^ = diag(l,l,e'^?/^,e-^^/^, 

gZA^/2T^ ^-iAl/2T^ ^iAl/2T^ ^-iAl/2T^ j^2 l^^ 

In a more general case, when 

Ao = A^T^ + A^qT^ , (2.14) 



the calculations of (|2.12p can be extended resulting in the 
modified form of the adjoint matrix; 

La = diag(l,l,e'^^/^,e-^^/^, 

^i{Al + V3Al)/2T ^-~i{Al+V3Al)/2T 
^'i{Al-/3Al)/2T^^-t{Al-V3Al)/2T\ ^ ^2.15) 



Since the rank of the SU(3) group is two, the Polyakov 
loop matrix in the adjoint representation La, can also be 
expressed in terms of two angler parameters 4>i and (f>2 
as 

La = diag(l,l,e*("^i-'^^) ,e-*('^l-'^=^ 



(2.16) 



Equation 112.121 1 can be generalized to any Nc- See e.g. [Till 
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Due to this change of variables, the partition function is with the coefficients C„ given by 
rewritten as 



\nZ = V 



(2^ 



Indet fl - L^e-IPl/^) + In M((/)i, (/^z) 



(2.17) 



with M being the Haar measure for a fixed volume V 
given by 



X sm 



2 

2(/)2 



sm 



which is normalized such that 







(2.18) 



(2.19) 







The first term in Eq. p.l7p yields the following ther- 
modynamics potential: 



2T 



(2^ 



trlnfl-LAe"^"/'^ 



(2.20) 



Mj the 



where in the quasi-gluon energy Eg = y 

effective gluon mass Mg is introduced from phenomeno- 
logical reasons To calculate the thermodynamic po- 
tential ()2.20p one still needs to perform the trace in a 
color space. We define gauge invariant quantities, nor- 
malized by dimensions of representations, as 



1 



N? - 1 



(2.21) 



where L^? is the Polyakov loop matrix in the fundamental 
representation, 



dia; 



iag (e'-^i 



g#2 g— 4(01+02] 



(2.22) 



and $4 is related with $ and $ via 





= 1, 


Ci 


= C7 




= Ce 




- Cs 




= 2 



1-9$*, 

1 - 27$<f> + 27 ' '^^^ 
-2 + 27l>$-8l($$) 



\ 2 



-1 + 9$$ - 27 + $^) + 81 ($$)^ 



(2.25) 

Thus, the n-body gluon contributions to the thermody- 
namic potential (|2.24p arc characterized solely by the 
Polyakov loop, i.e. the characters of the fundamental and 
the conjugate representations of the color SU(3) gauge 
group. 

The Haar measure (j2.18p is also expressed in terms of 
$ and $ as 



M(<^1,<^2) 

97r2 L 



, 6$$ -h 4 $3) -3($$)^ (2.26) 



A complete effective thermodynamic potential of glu- 
ons in the large volume limit is obtained from Eq. (|2.17p 
as follows: 



O = f7„ + 0$ + Co , 



(2.27) 



where fig is given by Eq. (|2.24p and the Haar measure 
contribution 



-floTln 



1 - 6$$ + 4 ($3 + $3) - 3 ($$)' 



(2.28) 



In (j2.28p we have neglected the normalization factor of 
the Haar measure which gives sub-leading contribution 
to thermodynamics. The ao and cq and/or gluon mass 
are free parameters which have to be fixed through cer- 
tain external conditions. They can be e.g. chosen to 
reproduce the equation of state of the SU(3) pure gauge 
theory obtained on the lattice through Monte Carlo cal- 
culations. 



(N^ - 1) $A = iV2$$ - 1 



(2.23) 



Carrying out the trace over colors and expressing it in 
terms of <I> and its conjugate one finally finds 



2T 



(2^ 



In 1 + ^ C„ 



-nEg/T 



(2.24) 



The SU(3) gluon plasma was studied in Q starting from the 
same thermodynamic potential as l|2.20| l. whereas the logarithm 
was expanded. 



3. ASYMPTOTIC EXPANSIONS OF THE 
POTENTIAL 

The potential for the SU(3) Yang-MiUs theory ^^TF!) 
obtained in the previous section provides an effective 
description of gluon thermodynamics. In particular, in 
asymptotically high temperatures, the fig should repro- 
duce the ideal gas limit. Indeed, taking the limit 
$, $ — > 1, corresponding to ~^ 0, one finds from 
dOl]), that 



17g($ = <I. = l) 



ler 



dP'p 
(2^ 



In 1 



-Eg/T^ 



(3.1) 
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Thus, the standard expression for a non-interacting gas 
of massive/massless gluons, is recovered. 

On the other hand, having in mind a quasi-particle 
approach, where gluons are considered as massive parti- 
cles with a temperature-dependent mass Mg = Mg{T), 
one can expand the logarithm in Eq. (|2.24|) . For a suf- 
ficiently large Mg{T)/T one approximates the logarithm 
by the first term of the expansion, resulting in the fol- 
lowing form of the potential; 



n„ ^ ^ 



J2^K,{nl3Mg)., 



(3.2) 



n=l 



where C„ are as in Eq. (|2.25p and Ki{x) is the Bessel 
function. The above can also be considered as a strong 
coupling expansion regarding the relation Mg(T) = 
g(T)T, where g{T) is an effective gauge coupling. 

The character expansion of the potential (|3.2p cor- 
responds to that obtained in the Polyakov loop mod- 
els on the lattice, derived using strong coupling tech- 
niques for the non-Abelian gauge group SU(3). Indeed, 
in the strong-coupling expansion the effective action to 
the next-to-leading order is obtained in terms of group 
characters as 



S, 



(SC) 
off 



AlO'5'lO + A2o5'20 + Alls'!! -I- A2lS'21 



(3.3) 



with products of characters Spq, specified by two integers 
p and q counting the numbers of fundamental and conju- 
gate representations, and couplings Xpq being real func- 
tions of temperature. One readily finds those couplings 
from Eq. p.2p . as well as the correspondence between 



Spq and C„ from Eqs. 



and (1123), as 



Ci,7 — Sio , C2,e — 5*21 , 

(^3,5 = 5*11 , C4 = 5*20 . 



(3.4) 



Taking only the contribution of a single-gluon distri- 
bution exp[— A/g/T] in the expansion Eq. (|3.2p yields the 
"minimal model" described by 



~T(T,Mg)<S><S>, 



(3.5) 



with the negative sign needed to get a first-order tran- 
sition as studied in [lO[. Here, an explicit form of the 
function J- relies on approximations used in evaluation 
of the momentum integration and parameterization of 
Mg{T). Assuming appropriate temperature dependence 
of so that the thermodynamic potential (|2.27p yields 
the phase transition with thermal expectation value of 
the Polyakov loop (<&) as the order parameter of Z(3) 
symmetry, the form widely used in the PNJL model |13l - 
UM is recovered. 



The action at this order includes one additional term, piVi with 
Vi = Y1s(\^f{^)\^ ~ 1) [loll . Since we deal with the temporal 
gluon field as a uniform background, Vi is not distinguished from 

•SlQ. 



In addition, the logarithm of the Haar measure part 
ri$ can also be expanded in powers of Z(3)-invariant op- 
erators. In this case, the effective gluon potential is found 
in the polynomial form 0, [l^ . Such form of ^l can, how- 
ever, be applied only to a weak first-order phase transi- 
tion. The polynomial form applied in the PNJL model 
was also shown to cause some problems in behaviors of 
charge fluctuations [l5| as well as with the phase struc- 
ture and symmetry properties of the potential at complex 
chemical potential |17l |. 



A HYBRID APPROACH 



The model described by Eqs. ((2:24)) . ((2:27| and ((2:28)) 
works fairly well when the thermal expectation ($) of the 
Polyakov loop is non- vanishing. However, at any temper- 
atures below Tc, where ($) = is dynamically favored, it 
causes some unphysical results on the equation of state. 

The model parameters ag and cq are flxed such that 
the model reproduces the value of Tc and the pressure 
at Tc calculated from SU(3) Lattice Gauge Theory. In 
such formulation, the model applied to the phase below 
Tc yields a positive pressure, however the entropy and 
energy densities turn out to be negative. In fact, keeping 
at low temperatures only the Ci term as a main contri- 
bution, one gets the potential 



Og(<I> = $ = 0) ~ 2r 



(2^ 



In I 



(4.1) 



which does not posses the correct sign in front of 
exTp[—Eg/T] expected from the Bose-Einstein statistics. 
One immediately flnds that the entropy and energy den- 
sities calculated from Eq. ()4.ip are negative. Therefore, 
the model cannot be naively applied to the phase below 
Tc- This problem appears not only with the complete po- 
tential ()2.24p , but also with its minimal form p.Sp which 
is frequently used in the literature. There, the equa- 
tion of state is entirely zero, at any temperature below 
Tc- This is clearly unphysical as there are color-singlet 
hadrons, glueballs, contributing to thermodynamics in a 
pure Yang-Mills theory and they must generate a non- 
vanishing pressure. 

The above aspects are in a striking contrast to the 
quark sector of the thermodynamics obtained in the pres- 
ence of a background gluon field Aq. There, the thermo- 
dynamic potential for quarks and anti-quarks with Nf 
flavors is obtained as. 



nq+q = -2NfT 



d^p 



(27r)3 



trln 



1 + Lpe 



-(E,-^i)/T 



+ (a« ^ -^^) 



(4.2) 



The trace over color indices in this case is easily per- 
formed and the potential is expressed by characters 
of fundamental $ and the conjugate $ representation 
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FIG. 1: Thermal expectation value of the Polyakov loop (left-hand figure) and the normalized pressure (right-hand figure) 
calculated in the hybrid model (|4.8p for massless Mg = and massive Mg = 850 MeV gluons. 



as 



n^+g = -2NfT 



X In 



(27r)3 



-E+ /T 



+ (m -m) 



(4.3) 



with E = Eq ^ n being the energy of a quark (anti- 
quark). 

In the limit of 3>, $ 0, which is expected at low tem- 
peratures, the contribution of one- and two-quark states 
is suppressed and only three-quark (baryonic) states, 
survives. This, on a qualitative level, 
is similar to confinement in QCD thermodynamics [l5| . 
One should however keep in mind that this model yields 
only colored quarks being statistically suppressed at low 
temperatures. On the other hand, unphysical thermo- 
dynamics that comes out below Tc from the gluon sec- 
tor (|2.27p . apparently indicates that gluons arc physically 
forbidden bellow Tc- We note that in the mean field ap- 
proximation higher representations of the Polyakov loop 
than the fundamental one are non- vanishing below Tc- 
This is an artifact of the approximation and in fact they 
do not condense when one evaluates their group aver- 
ages with the Haar measure. Such "hidden" physics of 
the higher representations can be embedded in the mean 
field approach when all the gluon energy distributions 
are expressed in terms of the fundamental Polyakov loop 
$. In this way all the colored gluons are suppressed and 
therefore the correct physics interpretation is recovered. 

This property is not affected by quarks. Let us con- 
sider massive gluons and quarks at zero chemical poten- 
tial. From Eqs. (|i?T|) and applied to T < Tc, the 
thermodynamic potential is approximated as 

'MA 2Nf 



y2 
2 



TT 



(4.4) 



Assuming that the glueball and nuclcon are made from 
two gluons and three quarks respectively and putting em- 
pirical numbers Mgiuobaii = 1.7 GeV m and M„ueieon = 
0.94 GeV, one finds Mg = 0.85 GeV and Mg = 0.31 
GeV. Given those numbers, Eq. (|4.4p is positive for ei- 
ther TV/ = 2 or 3. Consequently, the entropy density is 
negative at any temperature as found in a pure Yang- 
Mills case. Therefore, thermodynamics remains unphys- 
ical, unless additional terms responsible for the non- 
perturbative effects in confined phase are considered. 



4.1. Modeling glueballs as dilaton fields 

Below Tc, thermodynamics needs to be described in 
terms of physical degrees of freedom, i.e. glueballs. We 
introduce a glueball as a dilaton field x representing the 
gluon composite {AumA^'^), which is responsible for the 
QCD trace anomaly [201. The Lagrangian that we use is 
of a standard form given by 



1 



2 

4 VXo 



In ( ^ 



(4.5) 



where B is the bag constant and xo is a dimension- 
ful quantity. The two parameters B and xo can be 
fixed using the vacuum energy density S ~ jB = 0.6 
GeV fm""^ [2l[ and the vacimm glueball mass = 1.7 
GeV [l^ with the following definition: 



dx^ 



4:B 



(4.6) 



One finds that B = (0.368 GeV)" and xo = 0.16 GeV. 
With the Lagrangian (|4.5|) . the thermodynamic poten- 
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tial of effective glueball fields is found to be 




where a constant B/A \s added so tliat f2 = at zero 
temperature. 

4.2. The hybrid thermodynamic potential 

To avoid problems of unphysical equations of state in 
confined phase we adopt a hybrid approach which ac- 
counts for gluons and glueballs degrees of freedom by 
combining Eqs. (|2.27p and (|4.7p as follows: 

i7 = e(rc-T)o(x) + e(r-rc)o($). (4.8) 

The model parameters are constrained by requiring that 

• ri(fE') yields a first-order phase transition at Tc = 
270 MeV as found in SU(3) lattice calculations d, 

• Q.(x) and match at T^. 

When the gluon effective mass is assumed to be zero, one 
finds the following model parameters: 

Mg = : ($)t, 0.395 , oq (0.197 GeV)^ , 

Co = -(0.180 GeV)^ (4.9) 

One can also assume that the gluon becomes massive via 
non- vanishing gluon condensation (x) 7^ [2^. Requir- 
ing that a glueball is composed of two constituent gluons 
yields Mg = M^^jl. Since (x) little varies with tem- 
perature around f23 |. we treat Mg as a constant and 
choose Mg = 1.7GeV/2 = 0.85 GeV. The parameters in 
this case are found as 

Mg = 0.85 GeV : ($)t, = 0.439 , ao = (0.125 GeV)^ , 
Co = -(0.130 GeV)*. (4.10) 

A more general case, not considered in this paper, would 
include the temperature-dependent effective gluon mass 
which is fixed such that the present model quantifies ther- 
modynamics calculated on the lattice in the SU(3) gauge 
theory. 

4.3. Thermodynamics 

Thermodynamic properties of the hybrid model and 
its phase structure can be quantified directly from the 
potential (|4.8p . Fig. [T] shows the thermal expectation 



value of the Polyakov loop (4>) obtained from Eq. (|4.8p 
as the solution of the stationary condition, 951/9$ = 0. 
There is a trivial solution ($) = at any temperatures 
and it becomes degenerate with a non-trivial solution (<I>) 
at some Tc, indicating a first-order deconfinement tran- 
sition. 

The Polyakov loop expectation is weakly changing with 
Mg and approaches unity rather quickly as seen in Fig.[T] 
The temperature dependence of ($) just above Tc and 
its value at Tc are consistent with lattice results and can 
be still improved by introducing a thermal gluon mass 
Mg{T) as done e.g. in Q. There, the effective mass 
was parameterized as in the standard quasi-particle ap- 
proaches at high temperature, Mg(T) = g{T)T, with 
the effective running coupling g{T) Q. The lattice data 
on the renormalized ($) are known to exceed unity at 
T/Tc ^ 3 [2^. This property of lattice data, which is 
associated with uncertainties of the renormalization pro- 
cedure [1,[2^ , can never appear in effective Polyakov loop 
models where $ is the character of the fundamental rep- 
resentation and restricts the target space, so that (4>) is 
not allowed to go beyond unity. 

Fig. [T]-right shows the pressure calculated from the ef- 
fective gluon (|2.27p and from the effective glueball (|4.7p 
potential for massless and massive gluons. Although the 
presence of a constant cq in Eq. (|2.27p makes the pres- 
sure positive below Tc, as mentioned in the previous sec- 
tion, Eq. (|2.27p unavoidably leads to a negative entropy 
and energy densities. Consequently, in the hybrid model 
(|4.8p and for T < Tc, the pressure must be quantified 
by the glueball potential (|4.7p . The cusp at Tc in pres- 
sure implies a discontinuity in its temperature deriva- 
tive. The energy density and the interaction measure 
A ^ {E - 3P)/T^ are presented in Fig. H The energy 
density has a jump from glueballs to gluons thermody- 
namics at Tc, whereas the interaction measure exhibits a 
maximum just above Tc. 

Even though, the qualitative behaviors of E/T'^ and 
A follow general trends seen in lattice data, the EoS is 
apparently more sensitive to Mg than ($). The model 
calculations with massless gluons converge too quickly to 
asymptotic Stefan-Boltzmann limit. The EoS for massive 
gluons, on the other hand, exhibits a better agreement 
with lattice data. This clearly indicates the presence of 
some residual interactions above Tc which can be incor- 
porated into Mg{T). Thus, to quantify lattice results one 
would need to include the temperature-dependent gluon 
mass. In simplified quasi-particle models, where $ = 1 
for any T, the Mg{T) was shown to be strongly increas- 
ing when approaching Tc from above. This behavior, 
however, can be modified if the contribution of the back- 
ground gauge field is included in the quasiparticle model 
formulated in Eq. (|2.24p . 

For T < Tc the hybrid model includes a glueball as the 
relevant degree of freedom, similarly as seen in lattice cal- 
culations. However, in the present treatment, the model 
contains only the lowest-lying glueball, which might not 
be sufficient to quantify lattice thermodynamics in a con- 
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fined phase. One way out would be to deal with a gluon 
degeneracy factor as an additional parameter. 

5. CONCLUSIONS 

We have derived the thermodynamic potential in the 
SU(3) pure Yang-Mills theory in the presence of a uni- 
form gluon field within the background field method. We 
have shown that such effective gluon potential, which ac- 
counts for quantum statistics and reproduces an ideal gas 
limit at high temperatures, is formulated in terms of the 
Polyakov loop in the fundamental representation. 

The gluon distributions are found to be specified solely 
by the Polyakov loop and therefore there is one-to-one 
correspondence to the effective action in the strong- 
coupling expansion. We have shown that effective models 
of the Polyakov loop used so far to describe pure gauge 
theory thermodynamics appear as limiting cases of our 
result. 

Our main observation is that the effective Polyakov- 
loop potential can not be applied to the phase where the 
thermal expectation value of the Polyakov-loop vanishes. 
There, in confined phase, the equation of state is unphys- 
ical resulting in negative entropy and energy densities. 
This property of gluon potential is in remarkable con- 
trast to the description of "confinement" within a class 



of chiral models with Polyakov loops [isl . [T6j . There, 
colored quarks are suppressed only statistically at low 
temperatures. 

The gluonic model considered here indicates that col- 
ored gluons are forbidden below Tc as dynamical degrees 
of freedom. This feature is unchanged by the presence 
of quarks. To avoid problems of unphysical thermody- 
namics in confined phase we proposed a hybrid approach 
which matches at dcconfinemcnt critical temperature the 
effective model of gluons to the one of glueballs con- 
strained by the QCD trace anomaly. 

The approach developed in this work is open to fur- 
ther investigations of the SU(3) gluodynamics guided by 
available lattice results, as well as to a more realistic 
description of an effective QCD thermodynamics with 
quarks. 
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